Planar Spin Network Coherent States 
II. Matrix Elements 

Donald E. Neville 
Department of Physics 
^ I Temple University 

§ ■ Philadelphia 19122, Pa. 

(N' 



l> 



> 

in 



o 



X 



July 5, 2008 
Abstract 



^' This paper is the second of a series of three which construct coher- 

ent states for spin networks with planar symmetry. Paper 1 pro- 

bJO! poses a candidate set of coherent states. The present paper demon- 

strates explicitly that these states are approximate eigenvectors of 
the holonomy and momentum operators E} (as expected for coher- 
ent states), up to small correction terms. Those correction terms are 



cn . calculated. A complete subset of the overcomplete set of coherent 

^ ' states is constructed, and used to compute the approximate inverse 

of the volume operator. A theorem of Thiemann and Winkler is used 
Q ■ to calculate the matrix elements of the [volume, holonomy] commu- 

OO ! tator. In the classical limit this commutator takes the derivative of 

volume with respect to angular momentum L. 
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I Introduction 

This paper is the second of a series of three which construct coherent 
states for spin networks with planar symmetry. The first paper of 
this series proposed a candidate set of coherent states for planar 
gravity waves pT]. Sections HTl and UTTl of the_ present paper calculate 
the action of the densitized inverse triads E} and the spin one-half 
holonomy on these states. In each case, the leading contribution 
to each matrix element is derived; and qualitative estimates are 
given for the size of the small correction (SC) terms. Quantitative 
measures for the size of SC terms are worked out in the appendices, 
especially appendix [Dl The third paper in the series will discuss the 
behavior of the Hamiltonian in the classical limit |2] • 

The coherent states form an overcomplete set, and as such are not 
well suited as a basis for perturbation theory. Section HVl constructs 
a complete subset of the overcomplete set, and uses the subset to 
construct a series expansion for the inverse of the volume operator. 

The spin network Hamiltonian depends only on the E and the 
spin 1/2 holonomy, and matrix elements of these operators are 
worked out in sections [TTl and IIIII However, the Hamiltonian con- 
tains a commutator of volume with holonomy, suggested by Thie- 
mann to regularize a dangerous factor of inverse triad determinant 
in the field theoretic Hamiltonian [3J. This commutator is sensitive 
to the small correction terms (SC terms) and is difficult to compute 
using only the results of sections [TTl and lllli However, a result due 
to Thiemann and Winkler does the trick [4]. Accordingly (sections 
[TVl discusses the volume operator and its inverse, then) section |V] 
derives a formula for the commutator, using the Thiemann Winkler 
result. 

In appendix O I must restrict the magnitude p to be large, e^ :^ 1, 
say p greater than 5 or so, in order to obtain manageable approx- 
imations for the D(H) factors in the coherent state wavefunction. 
Whereas the matrix u determines peak values of the holonomy, 
therefore peak values of angle, the vector p determines peak value 
of the conjugate variable, angular momentum L. Specifically, from 
eq. fl85|) . the magnitude p determines the peak value of the magni- 
tude L of the angular momentum. 

<L + l/2 >=p/t. 

Since we are focusing on large < L >, the restriction to large p is 



not a problem. 

The following is a summary of the most important formulas from 
paper 1. The proposed coherent states are sums of 0(3) spherical 
harmonics. 

|M,p) = iV^((2L + l)/47r) 

■ exp[-tL(L + l)/2]D^''\h\MD^'^\g\)M^. (1) 

The rotation matrix D(h) is a spherical harmonic, apart from nor- 
malization. 

v/47r/(2L + l)YLM(^/2,</.-7r/2) = D(L)(h)oM(-0+vr/2, 0/2, </.-7r/2). 

(2) 
The matrix g, eq. ([1]), is in SL(2,C); every such matrix can be de- 
composed into a product of a Hermitean matrix H times a unitary 
matrix u. 

g'-^^ = exp[S ■ p]u'^^^ 

:= i7(^)«(^) (3) 

All three matrices (h, H, and u) have their axis of rotation in the 
XY plane. 

(4) 
(5) 



/.(^) 


= exp[im ■ S'6'/2]; 


m 


= (cos 0, sin 0,0). 


w(^) = 


u{-p + 7r/2,a,P-n/2) 


:= 


exp[in ■ S a/2]] 


h = 


(cos/?,sin/3,0). 


((L) = 


exp[S' ■ p]] 


P = 


p [cos(/? + fi), sin(/5 + /i), . 



(6) 



(7) 

From the qualitative discussion given in paper 1, the matrix u should 
determine the peak value of the holonomy h, while p should deter- 
mine the peak angular momentum. 

The Hamiltonian is expected to contain both densitized triads E 
and holonomies /i^^/^^ in the fundamental representation of SU(2). 



However, it is desirable to replace the SU(2) representations in 
the Hamiltonian by 0(3) harmonics, and construct the coherent 
states as a series of 0(3) harmonics, because the action of the E 
on 0(3) harmonics is especially simple; see for example section [III 
of the present paper. In particular, matrix elements of the SU(2) 
holonomies in the Hamiltonian, hmn , should be replaced by 0(3) 
spherical harmonics, by using the formulas 

Vih)M = D^'\h)oM. (8) 

From the last line of eq. ([8]) the components of the unit vector V{h) 
form an L = 1 spherical harmonic of 0(3). In section UTTl on matrix 
elements of the spin 1/2 holonomy, I replace the holonomy by V{h), 
in effect replacing SU(2) by 0(3). This replacement makes it easier 
to compute the action of the holonomy on the 0(3) coherent states. 
Eq. (Q can also represent the angular wavefunction of a particle 
moving in a central potential. In this interpretation, V{h) is the 
unit radius vector of the fictitious particle, and V{u) is the peak 
value of the radius vector. 

II Action of the E 

This section computes the action of an E operator on the coherent 
state. The main results of this section are the following: the coherent 
states are approximate eigenvectors of the E , 

(7fi;/2)^^E^ I u, p) =< L > (ua coSyU — h x V)Asin/i) | u, p) + SC. 

< L > is the average, or peak value of the angular momentum, n is 
the axis of rotation for u, eq. ([H]), while /i is the angle between fi and 
P, eq. ([7]). V is the vector V{h), eq. ([8]), evaluated at the value h = 
u. (The next section on the holonomy will show that V is the peak 
value of V{h), hence u is the peak value of h.) The small correction 
terms SC are investigated in appendix [D] and shown to be down by 
factors of order 1/a/L. 



I now derive eq. ([9]). By construction, the 0(3) D(h) matrices in 
the coherent state transform simply under the action of an E : 

{ik/2)-'W^ I w,p} = iv5^((2L + l)/47r) exp[-tL(L + l)/2] 

L,M 

■D^''\h)oN{L,N I Sa I L,M)D^'^\g^)Mo 
= N J2 ((2^ + l)/4vr) exp[-tL{L + l)/2] 

L,Af,R 

= iV ^((2L + l)/47r) exp[-tL(L + l)/2] 

•[D(^)(/i)S^D(^)(«t)D(^)(/7)]oo (10) 

On the second hne I use the g = Hu decomposition, eq. ([3]), plus 
D{g^) = D{u])D{H). On the last line I have used a matrix notation 
to hide some indices. 

I now implement the general procedure outlined in section IV 
of paper 1. This procedure involves three steps. Produce a factor 
of D^^\u), because the rows of this matrix provide a convenient 
basis of unit vectors; rewrite the summand until it looks as much 
as possible hke the summand in the original coherent state, because 
we wish to prove that \ u,p) is an approximate eigenfunction of E ; 
simplify the summand by power series expanding and integrating 
term by term, because D(H) is Gaussian, and higher powers in the 
expansion will be suppressed. 

First produce a factor of D^^^ (u). The spin generator S in eq. flTOj) 
is essentially a Clebsch-Gordan coefficient, and from the rotation 
properties of these coefficients, appendix [XJ eq. (jHSD, 

SADW(ut) = D(L)(ut)SBD«(u)BA. 
Eq. fITOl) becomes 



(7ft:/2)-iE^ I M, p) = nY^H^L + l)/47r) exp[-tL{L + l)/2] 



L,AI 

■[D(^)(M)SBD(^)(i7)]oo/^(^Hw)BA 
= A^ '^{{21 + l)/47r) exp[-tL{L + l)/2] 

L,M 

•l^(''nM)oM/^(''n^)Mo[(I^,M I Sb I L,M-5)D(^)(/7)m-b,o//^^''H^)mo] 

■I^('n«)BA. (11) 

The final two lines exhibit the desired factor of D^^\u), as well as 
factors of D^^\hu\)oM D''^\H) mo which make the sum in eq. (ITT|1 
look as much as possible like the original coherent state. 

I must now evaluate and expand the square bracket, second line 
from the end of eq. (TTT]) . The ratio of D(H) matrices follows from 
eq. (176|) of appendix [C| vahd for e^ ^ 1. 



D(l)(H)m^i,o = D(l)(H)m,o[1/(pi - ip2)]^V(L±M)/(LTM+l). 

(12) 
I have specialized to p^ = 0. For this value of pa I can parameterize 
p as at eq. (I7l), so that pi — ip2 = exp[—i{P + /i)]. Also, the matrix 
element of Sq is just M; while Sb, B = ±1, is given by 



(L, M I S±i I L, M T 1) = Tv^(L T M + 1)(L ± M)/2 (13) 

Eq. (TTT!) becomes 



{^k/2)-'W;,\u,p) = N J2 ((2L+l)/47r)exp[-tL(L + l)/2] 

L,M,B 



C[M,B = 0] = M- 
C[M,S = +1] = -{L + M)/V2- 
C[M,B = -l] = {L-M)/V2. (14) 

I replace the D^^\u)ba by the geometrically more transparent 
quantities n, V, and n x V. The row D^^\u)oa is just the vector 
V{h) introduced at eq. (jHD, with h replaced by u. 



D«(u)oa = V(u)a. 

The rows D^^\u)±i^a may be replaced by linear combinations of the 
unit vectors h and h x V, using eqs. (1^ and (17D1) in appendix [Bl 
Eq. (HM becomes 

(7«:/2)-i% I w,p} = iV5^((2L + l)/47r) exp[-tL(L + l)/2] 

L,M 

■D^'^\hu\),MD^'^\H)Mo 
■ {(Lcos/i + iMsin/i)n^ + (— Lsin/i + iM cos fi){n x V^)^ 

+ My(M)A}. (15) 

The third and final step of the general procedure invokes the 
Gaussian nature of the factor D(H). From appendix [Cl 

D(^)(i7)Moexp[-tL(L + l)/2] = 

^/tJ^ exp[-t{L + 1/2 - p/tf/2] 

Vl/[(L + l/2)7r] exp{-MV[2(L + 1/2)]} f{p, t). 

This suggests expanding the factors of L and M in the curly bracket 
around their means, < L >= p/t -1/2 and < M > = 0. Equiva- 
lently, in eq. flTSl) set 

L =< L > +AL; M = AM, 

and keep out to first order in AL and AM. The dominant terms in 
the expansion come from the terms in the square bracket involving 
< L >, 

< L > (ii cos /i — h X V) A sin /i) . 

This expression is a constant and may be taken out of the sum, 
leaving behind just the original coherent state. We get the desired 
result eq. (^ for the leading term. 

The remaining terms in the expansion are AX terms, X = L or 
M. These small correction (SC) terms are studied in appendix [Dl 
where it is shown they are down by factors of order ax/ < L >, ax 
the standard deviation for the Gaussian distribution of X. 



Ill Matrix Elements of the Holonomy 

This section computes the action of a holonomy /i^^/^^ on a coherent 
state. The main results are the following. When h acts on a coherent 
state, all of the states which result differ from the original coherent 
state. In particular, the leading terms are not identically equal to the 
original ket. To expand the action of the holonomy, it is necessary 
to introduce new kets 



|L±1) = iV(L ± 1) ^((2L + l)/47r) exp[-tL(L + l)/2] 

L,M 

■D^'^^'\huf)oMD^''\H)Mo. (16) 

The first part of this section shows that these kets are Gaussian 
distributed in L, but with peak values of L shifted by one unit, from 
< L > to < L > ±1. (Strictly speaking the ket notation should 
be < L > ±1 rather than L ± 1.) The second part of this section 
expresses the action of the holonomy in terms of these states and 
shows that 



D«(h)oA I u,p) = D«(u)oa(| L + 1)+ I L - l))/2 + SC. (17) 

From eq. ([8]), the 0(3) spherical harmonics Z^^^-*(/i)oa are linear com- 
binations of /i^/^ matrix elements. Eq. fITTI) states that D^^\u)oa is 
the peak value of D^^\h)oA', consequently u is the peak value of h. 
The linear combination (|L + 1)+|L — l))/2is approximately the 
original ket \ u,p). 

Section III of paper 1 describes a electron-in-a-central-force anal- 
ogy which predicts that the expectation values of holonomy and E 
should be orthogonal. 

= ^<DW(h)oA><E^>. (18) 

A 

From eq. ([8]), D^^\u)oa in eq. (TTTj) is just Va. Therefore the con- 
straint given in eq. (TTSi) is obeyed, since < E >, eq. ([9]), contains 
only vectors perpendicular to V. 



A the kets | L ± 1) 

The states defined at eq. (fT6!) have peak angular momentum shifted 
by one unit. Proof: in eq. (TTB]) . the D(H) factor is easier to approxi- 
mate, while D{hu]) is harder. Therefore relabel L±l = L, then drop 
tildes, in order to make the D{hu\) factor equal to the correspond- 
ing factor in the original coherent state | M,p): under the relabeling, 
D^^^^\hu\)-^D^^\hu\)^D^^\hu\). The relabeling changes the 
exponential exp[— tL(L + l)/2] significantly; and it produces negli- 
gible changes in (2L -|- l)/47r = (2L + l)/47r -|- order(l/L) The new 
D(H) factor is D^^\H)^D^^'^^\H)^D^^^^\H). This matrix dif- 
fers from the corresponding matrix in | u^p), but is easy to relate 
to that matrix. I group together the two factors in eq. (TT6|) that 
change significantly, and use eq. (185|) . 

exp[-t(L ^ 1)(L T 1 + l)/2]D(^^i)(iJ)Mo 
^ exp{-t[L T 1 + 1/2] V2 + t/8} exp[-iM/3](exp(p/2)/2)2^^2 

{2Lt2)\ 

^{L T 1)!^(^ T 1 + M)\{L T 1 - M)! 

^ exp{-t[L ^ 1 + 1/2 - p/tf/2 + p^/2t - p/2} 

exp{-MV[2(L T 1 + l/2)]}/v/7r(L ^ 1 + 1/2) 

= exp{-t[L T 1 + 1/2 - p/t] V2 + P^/2t - p/2} 

exp{-MV[2(L + l/2)]}/v/7r(L + 1/2). (19) 

The factorials were replaced by a Gaussian using eq. (ISTl) . Now 
compare the Gaussian in eq. (TT^ to the Gaussian in the original 
coherent state, eq. fl85j) . The two Gaussians in L are the same 
except for a shift in peak value. 

new < L + 1/2 >= p/t ± 1 = original < L + 1/2 > ±1. (20) 
D 

B Action of the holonomy 

Presumably the Hamiltonian will be a function of the matrix h^^^"^^; 
but for calculations it is more convenient to use the linear combina- 
tions D'^^\h), eq. (jH]), rather than the h^^/'^\ The coherent state is 



an expansion in spherical harmonics D^^\h)^ and using the matrix 
D^^\h) allows one to invoke identities for simplifying the product 
of two D's. 



D(i) {h)oA I M, p) = iV ^((2L + l)/47r) exp[-tL(L + l)/2] 

L,M 

•5^D(^±i)(/i)oAf(I^±l,0|L,0;l,0) 

■ (L±1,M I L,M-A;1,A) 

■D^''\u])M-A,ND^''\H)m (21) 

On the right I have rewritten D(g\) = D{u\)D(H), from eq. ([3]); 
and I have used the formula for combining two D(h) matrices. 



L±l 

• (L ± 1, I L, 0; 1, 0)(L ± 1, M I L, M - A- 1, A). 



This formula is a special case of the general relation for coupling 
two D's, eq. (!64|l . Note the value L is excluded from the sum on the 
right. (L', I L, 0; 1, 0) vanishes for L' = L. Parity conservation also 
predicts no contribution from L' = L. 

I now implement the same three step procedure as used in the 
previous section for E matrix elements. I introduce the basis vectors 
D^^'iu) by invoking eq. (1641) . 



(L ± 1, M I L, M - A; 1, A)D^''\u])m-a,n 

= D^'^^'\u])m,n+b{L±1,N + B \ L,N;1,B)D^'\u)ba. 
I insert this into eq. (I2T!) . and relabel N —^ N — B. 



10 



D^'\h)oA\u,p) 

= AT ^((2L + l)/47r) exp[-tL{L + l)/2] 



L,N 



■ [(L±1,0 I L,0;l,0)(L±l,iV | L,N -B;l,B) 

■ D'^''\H)n.b,o/D{H)n,o]D^'\u)ba. (22) 

I have multiplied and divided by D{H)no, so that the second and 
third lines resemble the original coherent state, except for the change 
L ^ L ± 1. 

I must now evaluate and expand the square bracket on the last 
two lines. I use eqs. flTTj) and fl76l) with p3 = to evaluate the ratio 

D(H)n_b,o/D(H)no. 

I relabel N — ;> M, and use the results of appendix [Bito replace D^^^ (u) 
by n, V, and fi x V. 

= N^{2L + l)/(47r) exp[-tL(L + l)/2] 



L,M 

I MO 



L±l 

■ {c{L ± 1, V)V{u)a + c{L ± 1, n)nA 

+ c{L±l,hxV){hxV)A). (23) 
The coefficients c are 
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c{L + l,V) = (L + 1,0 I L,0;1,0)(L + 1,M | L,M;1,0); 
c{L + l,n) = (L + 1,0 I L,0;1,0) 

■ J2iiTl/V2){L + 1, M I L, M T 1; 1, ±1) 



• exp{±ifx)^/{L ± M)/{L tM + 1)); 

c{L + l,nxV)=i{L + 1,0 mO;l,0) (24) 

•5^((-i/V2)(l + i,m|l,Mt1;1,±i) 
± 

• exp(±i/i) v/(i^ ± M)/{L tM + 1)). (25) 

For L+1 -^ L-1, replace (L + 1, M | and(L + 1, | by (L - 1, M | 
and(L — 1,0 |. After some work with a table of 3 J symbols, one 
finds 

c{L + i,v) = y/{L + iy-My{2L + iy, 

c{L-l,V) = VL2-M2/(2L + 1); 

c{L + l,n) = -[y/{L + 1 + M)/{L + 1 - M)(L + M)/(2(2L + 1)) 
-(M^ -M)]cos/i 

-2[v/(i: + 1 + M)/(L + 1 - M)(L + M)/(2(2L + 1)) 
+ (M^ -M)]sin/i; 
c(L-l,n) = z[VL2-MV(2L + l)]sin/i; 
c(L + l,nx1/) = -i[y^{L + 1 + M)/{L + 1 - M)(L + M)/(2(2L + 1)) 
+ (M-> -M)]cos/i 

+ [V(^ + 1 + M)/(L + 1 - M)(L + M)/(2(2L + 1)) 
-(M^ -M)]sin/i; 
c(L-l,nx\>) = i[VL2-M2/(2L + l)]cos/i. (26) 

Continuing with the general procedure, I expand the L and M 
dependence of the c's. I replace 

L->< L > +AL; M^AM. 

and keep terms in eq. ( l26l) out to linear in AX, X = L or M. 

12 



c{L + l,V) ^ <(L + 1)/(2L + 1) >; 

c(L-l,\/) ^ <L/(2L + 1)>; 

c{L + l,n) = -(M/ < L >)coSyU-2sin/i < L/(2L + 1) >; 

c(L-l,?T,) = zsin/i < L/(2L + 1) >; 

c{L + l,nxV) ^ -2 < L/(2L + 1) > cos/i + (M/ < L >)sin/i; 

c{L-l,nxV) ^ +i < L/(2L + 1) > cos/i. (27) 

Factors of L inside brackets <> should be interpreted as < L >; 
terms of order AX/ < L > have been kept, but not terms of order 
AX/ < L >^ or {AX/ < L >)^. After inserting this expansion into 
eq. (123|) . that equation becomes 



D^'\h)oA\u,p) = Ar^((2L + l)/47r)exp[-tL(L + l)/2] 

L,M 

■{(l/2)[D(^+i)(M) + D^'^-'\hu])]oMVA 

•[— cos fifiA + sin/i(n x V)a] 
-^(l/2)[Z}(^+i)(M) - D^''-'\hu])]oM 
■[siiafihA + cos /i(n x V^)a] 
+order((M/L)2, 1/L)}D^'^\H)mo. (2^ 

I anticipate that 



M/ < L >= order ctm/ < L >= order l/v< L >. 



Therefore I am keeping order 1/a/< L > but neglecting order {M/ L^ 
= order 1/L. Therefore when going from eq. (1271) to eq. (125]) it is 
legitimate to replace < L/{2L + 1) > by 1/2, etc.. 

Now separate eq. (I28p into a leading term plus small correc- 
tions. Eq. (1281) involves the linear combinations {1 / 2)[D'^^^^\hu\) ± 
D'^^^^\hu\)\, therefore linear combinations (1/2)(| L + l)± | L — 1)) 
of the kets investigated earlier in this section. The linear combina- 
tion with the upper sign is a sum of two Gaussians in L. The two 
Gaussians have slightly different means, but large standard devi- 
ations a^ = l/y/t, from eq. ( fT9l) . Therefore the two Gaussians 
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strongly overlap with each other and with the Gaussian of the orig- 
inal coherent state \u,p). Therefore, inside a sum over L, 

{l/2)[D^^+^\hu\) + D^^-^\hu\)]D^^\H) = D^^\hu\)D^^\H)- 
(1/2)(|L + 1)+|L-1)) = \u,p). (29) 

The linear combination with the upper sign will yield the dominant 
term. 

The linear combination with the lower sign, on the other hand, 
will be the difference of two very similar Gaussians, therefore will 
be very small. A detailed calculation in appendix [D] shows it to be 
a small correction. 

Also, the ft terms involve a factor of M/< L >, and the discussion 
of appendix [D] shows these terms are down by l/vL, as expected 
because D^^'{H) is Gaussian in M. Therefore only the Va term is 
dominant. From eq. ([H]) Va is also D^^\u)oa', and we get precisely 
eq. (HZD. 

In the leading, V term of eq. ( !28l) I have not replaced the linear 
combination (1/2)(D(^+^)(M) + D^^-^\hu])) by D^^^hu]), even 
though (1/2)(| L + 1)+ I L — 1)) overlaps strongly with | M,p). 
(Similarly, I have not replaced (1/2)(| L + 1)+ \ L — 1)) hj \ u,p) in 
eq. (IT7I) .) The linear combination (1/2)(| L + 1)+ | L — 1)) has the 
same E eigenvalues as | u,'p). However, (1/2)(| L + 1)+ | L — 1)) 
does not behave like | u, p) under the action of a holonomy: 

D«(%a(1/2)(|L + 1)+|L-1)) 

^(l/4)(|L + 2) + 2|L)+|L-2)), 

where | L) is the original state \ u,p). I am not yet sure whether 
these changes in peak L have any significance, and for now, I prefer 
to keep them explicit. Note if the Hamiltonian contains a product 
of n holonomies, then the original coherent state will be changed 
into a linear combination of coherent states with peak L shifted by 
as much as ±n. 

So far I have carried out all calculations assuming the axis of 
rotation for the matrix u and the vector p both lie in the XY plane. 
What happens if I replace (u, p) by parameters {u,p) not obey- 
ing this constraint? I state the following theorems without proof. 
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The new and old matrices H and H are connected by a similarity 
transformation. 

H = vfHv. 

The matrix v is uniquely determined by p and p, up to an overall 
sign: 

v = exp[io-- (p X p)x/2]. 

X is the angle between p and p. The new expectation value of D^"^' {h) 
will be 

<D«(h)oA>=D«(vii)oA, 

which means u cannot be chosen arbitrarily; vu must have its axis 
in the XY plane. The new expectation value of E will have the same 
value as before, except replace u by vm. 

All these theorems are relatively easy to prove: at an early step, 
use the similarity transformation to replace the poorly understood 
H by H; then use the extensive results of appendix [C] to approximate 
H. Note that the input unitary matrix u no longer gives the peak 
value of the holonomy h. 

Moving the parameters out of the XY plane seems to add nothing 
but complexity. If we stick with (u, p), we already have enough 
parameters to tune the expectation values of both holonomy and 
triad. For p^ = 0, 



<D«(h)oA>=D«(u; 



OAi 



therefore we can adjust the holonomy to any value by adjusting u. 
There are only two adjustable paraineters left, (pi,P2)- However, 
the remaining expectation value, < E^ >, must lie in a plane per- 
pendicular to D^^\u)qa- Therefore only two parameters are needed 
to tune its value. 



IV V sqared, V, and V inverse 

This paper is built on two approximations. The first approximation, 
an expansion in inverse powers of e^, is used in appendix Oto obtain 
a manageable expression for D(H). 

The second approximation is an expansion in powers of l/-\/< L >: 
appendix [Dl shows that small correction (SC) terms are down by fac- 



tors of 1/V< L >. 
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This latter expansion is not a traditional perturbative expan- 
sion: there is no Hamiltonian with an unperturbed part and known 
solutions, plus a small perturbing potential. However, there is a 
small parameter, 1/a/< L >; and perhaps more importantly, there 
is a complete set of orthonormal states. In this section I use this 
framework to compute the approximate square root and approxi- 
mate inverse of the volume operator. 

The next subsection needs a fact established in appendix [Dl the 
SC states are proportional to moments of the Gaussians in D(H). 
I. e. the original coherent state contains D(H), whereas the SC states 
contain D(H) times powers of M or (L - < L >). This is perhaps 
plausible because of the way the SC states are generated: by power- 
series expanding a summand in powers of M and (L - < L >). 

A SC States as a Complete Subset 

The expansion 

E I u, p) =< E > I u, p) + SC 

resembles a perturbation expansion in certain respects. It is an 
expansion in a small parameter, l/-\/< L >, since SC terms are down 
by l/\/ < L >. Also, the states contained in the SC terms, plus the 
original state | M,p}, form a complete subset of the set of coherent 
states. The set of coherent states 

{| u,p),Vu,p} 

is overcomplete, therefore not suitable for perturbation theory, which 
requires a complete set. The complete subset is therefore a more nat- 
ural basis for calculations than the original set of coherent states. 

The complete subset is constructed by repeated application of the 
E operators to the original coherent state and the states contained 
in the SC terms. Schematically, 

E |u,p) = <E >\u,p) + ^hix\ml{X))- 

X 

E |ml(X)) = <E >|ml(X)) 

+b I m2(M)) +c\u,p) 

+d| ml(M),ml(L)), (30) 
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etc.. The ket | ml{X)) is the original ket, | u,p), except that the 
Gaussian D^^\H)mo factors are replaced by 

(X-<X>)D(l)(H)mo, 

X = L or M; i. e. the ml stands for "first moment" of the Gaussians. 
Similarly, | m2{X)) is proportional to the second moment, (X— < 
X >y, etc.. (There is a similar series of SC states associated with 
the holonomy operator. In this section I consider only the series 
generated by the_E , because the operators discussed in this section 
depend only on E .) The bix and b, c, d factors are down by 1/vT; 
the terms containing these coefficients are small corrections. 

The presence of a "c" term in eq. (!30!) is perhaps surprising: the 
E operator has the ability to remove one power of M as well as add 
one power. For a full discussion, see eq. (JUj) in appendix [Dl 

Each state in the set 

{| mp(X), mq(Y)),V integers p,q} (31) 

can be shown to be a linear combination of states in the original 
coherent basis. (This is clear from the overcomplete nature of the 
original basis; and see appendix [D] for an example.) Therefore the 
series forms a subset of the original basis. 

This subset is complete and can be made orthogonal. On going 
from the original state to the SC states, one replaces the D(H) factor 
by D(H) times a power of X— < X >. Since the original factor D(H) 
is Gaussian in X, one is replacing the Gaussian by Gaussian times 
a Hermite polynomial, i.e. by a wavefunction for a simple harmonic 
oscillator (actually, two wavef unctions, one depending on M, and one 
depending on L, since D(H) is a Gaussian in M times a Gaussian in 
L). When the dot product of two SC states is computed, the factors 
depending on D(h) and D(u) always drop out. (Examples of this 
behavior are given in appendix [Dl) What remains is a sum over L 
and M, times products of simple harmonic oscillator wavefunctions 
in L and M. The sums over L and M can be turned into integrals, 
and the orthonormality follows. 

The complete subset is a natural basis to use when computing the 
square root of (^^V)^ and the inverse of *^^V . In particular, *^^V has a 
kernel, so that its inverse can only be approximate; the perturbation 
approach is useful in determining when the approximation breaks 
down. 
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B ((2)V)2 

This section calculates the action of the volume squared operator 
on the coherent states. This is the first point where I must take into 
account the fact that there is one holonomy for the x direction and 
another for the y direction. 

I U, p) ^1 Ux,Px;Uy,Py) =1 Ux,Px) I %,Py)- 

The volume operator contains both E;^, which acts only on the first 
ket, and E^ , which acts only on the second ket. 

The E matrix is block-diagonal because of gauge choices j5] . The 
volume operator *^^V therefore simplifies to 

:= E| (2)E _ (32) 

The holonomies along the z direction are eigenstates of E|. Since 
E| is diagonal in the basis chosen, the non-trivial part of the volume 
operator is 

((2V)2 = (2)E , (33) 

which is the determinant of the transverse components of E . In 
what follows I will refer to this (or its square root) as the "volume" 
operator, for short, even though "^^Vhas the dimensions of an area. 
(I cannot call '•^•V, e.g. transverse area; E^ is the transverse area.) 
The action of the E on the coherent states is given at eq. 



(7fi;/2) ^((^V)^ I U^,Px) I Uy,Py) = Aoo I Ux,Px) I Uy,Py) 

+ y^^ Aqix I Ux,Px) I ml{Xy)) 

Xy 

+ ^AlOX I ml{Xx)) I Uy,Py) 

Xx 



+ order l/^LlLl, p+q = 2. (34) 

The lowest order SC states are ml states (first moment states) such 
as I ml{Xy)) with Xy = My or Ly. From eq. (l88l) . the leading term 
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is order L^Ly, 

Aoo = {-f^i/2)-^<El><El>ezAB 

~ < Lx >< Ly > X order unity. (35) 

From eqs. (JHID and (11031) the non-leading terms are suppressed by 
factors of 1/a/L. 



A„,nx = order < L.Ly > (W < L, >r{\l/ < Ly >)°. (36) 



The action of (^^V)^ on other members of the complete set is similar, 
for example 

(7«:/2)-2((2V)2 I ml(M,)) | Uy,Py) = Aqo | ml(M,)) | Uy,Py) 

+ J^^oiMX I ml{M,)) I ml{Xy)) 
Xy 

+ BlQM I m2(M)) I Uy,Py) 

+ C'lOM I «x,Px) I Uy.Py)- (37) 

The -Bmn and Cmn are suppressed by the same factors as the Amn, 
eq. (l36ll . Eq. (!37j) implies that the other members of the complete set 
are also approximate eigenfunctions of (^^V)^with the same eigen- 
value (same leading coefficient Aqq). 

The above results suggest that ^^V , the square root of the volume 
operator, is also a double power series in 1/a/Z^, i = x,y . The series 
for ((^V)^ , eqs. f l34|) and fl37j) have L^Ly as highest power, which 
suggests that the '^^V series has highest power ^L^Ly. In the next 
subsection I propose an ansatz for *^^V having this form. 

C The Operator (^V 

The ansatz for *^^Vis 
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(7fi;/2) ^(^V I Ux,Px) I Uy,Py) = aoo I Ma;, Pa;) | Uy,Py) 

+ X^aoix I Ux,Px) I m^Xy)) 
Xy 

+ ^aiox I ml(Xa;)) I Mj/,Py) H ; 



Xx 



(7/t/2)~^(2V I ml(Xa;)) I My,Py) = aoo I ml{X^)) \ Uy,Py) 

+C10X I Ua:,Px) \Uy,Py) -\ ; 

(7/t/2)"^(2V I Ux,Px) I rnl{Xy) = aoo I u^,Px) \ ml{Xy)) 

+ COIX I U^,Px) \Uy,Py) -\ , (3^ 



where 



amn, c„,„ = ordervd:;(A/l7L;)"X\/l/Ly)"- (39) 



I wish to compute the a's, the coefficients giving the action of E 
on the original coherent state, to zeroth and ffist order. To do this 
I will need to know only the zeroth order term in the action of E 
on the higher moments; however, I will include the ffist order CmnX 
terms to make the point that the diagonal elements 

(Ux,Px I (%,Py I ^^V I Ux,Px) I %,Py) 

have no ffist order corrections. 
To determine the a's, I set 



(('V)2 I Ux,Px) I %,p,) = (^V ■ (2V I Ux,Px) I %,p,). (40) 

The left-hand side will be an expansion linear in A's; the right- 
hand side will be an expansion quadratic in a's. Both sides will 
be expansions in the orthonormal basis consisting of the original 
I Ux,Px) I Uy,Py), plus states with at most one | Ui,pi) replaced 
by a I ml{X)). (States with one m2 ket, or both original kets 
replaced by ml kets, do not contribute if we are keeping only zeroth 
plus ffist order.) As in standard perturbation theory, one can use 
orthonormality to project out the coefficient of each state. This 
yields a set of equations relating A's to expressions quadratic in a's. 
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The quadratic equations can be solved for the a's; phase ambiguities 
are resolved by demanding that ^^Vhave positive eigenvalues. 

Now write out both sides of eq. PUI) . The left-hand side of this 
equation was written out at eq. (1M|) . From eq. (155]) . the right-hand 
side is 



rhs = (^V[aoo | U^,Px) \ Uy,Py) 

+ '^aoix I u^,Px) I ml^Xy)) 

Xy 

+ ^aiox I ml(X^)) \uy,py) ^ ] 

Xx 

= flOO [OOO I Ux,Px) I Uy,Py) 

+ J^ooix I Ux,Px) I mli^Xy)) 

Xy 

+ ^a,iox I ml{X^)) \uy,py) ^ ] 

Xx 

+ y^Qoix[aoo I u^.Px) I m^{Xy)) 
Xy 

+C0IX I Ux.Px) \Uy,Py) ^ ] 

+ ^aiox[aoo I ml{X^)) \ Uy,Py) 

Xx 

+ CIOX I Ux,Px) \Uy,Py) -\ ] 

+ ■■■. (41) 



Set this expression equal to eq. (l34l) and invoke orthonormality to 
equate corresponding coefficients. 

(±1)^00 = Qqo + '^[QqixCqix + QioxCiox]; 

X 

{±l)Aoix = 2aooa.oix; 

i±l)Awx = 2aooaiox- (42) 

The ±1 = sign(y4oo) is inserted to insure that the quantity (±l)y4oo 
is positive. The matrix needed for construction of the Hamiltonian 
is I '•^V I, not *^^V; i.e. (*^^V)^may have eigenvalues of either sign, 
but the eigenvalues of *^^Vmust all be positive (or zero). Since the 
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diagonal, Aqq element of the ((^V)^ operator is largest, the sign of 
the eigenvalue will be determined by the sign of Aqq. 
The solutions to eq. ( H2l) are 

aoix = ±^oix/(2a/±A)o); 

aiox = ±Aiox/{2y^±Aoo). (43) 



Note the sum on the first line of eq. (H2l) is second order and may 
be neglected; therefore all three lines of eq. (1431) are already correct 
to first order. The first order corrections to Ooo vanish. 

D The Inverse of (^V 

The matrix "^^Vhas a kernel, and cannot possibly have an exact 
inverse. However, if '■^Vacts on a coherent state which is an ap- 
proximate eigenvector of *^^V , with very large eigenvalue, then one 
can neglect the states in the kernel, because *^^V connects the orig- 
inal state only to other (coherent and mn(X)) states having large 
eigenvalue: note the structure of eq. fl38l) . 

In this situation it is possible to calculate an approximate inverse 
for *^^V . As in the calculation of '•^V , I work only to zero and first 
order in 1/vL. 

The matrix elements of "^^Vwere assumed to be order ^yL^Ly, 
times a power series in powers of l/y/JTi. This suggests a corre- 
sponding ansatz for matrix elements of C-^V)^^: they are order 



1/ ^jLr^Ly^ times, again, a power series in powers of 1/ ^/Ll. 
Also, the approximate inverse must obey 

5(ix.jx)(5(iy,jy) = Snx,ny(jx,jy | ((^V)^^ | nx,ny)(nx,ny | (^V | ix,iy). 

(44) 
The indices (ix, iy), (jx, jy), (nx, ny) label states in the complete 
set, eq. (1311) . 

As a first step I keep only the leading, zeroth order, contributions 
to the matrix "^^V . In zeroth order *^^V is diagonal: 

(nx,ny | (^V | ix, iy) = (5(nx, ix)(5(ny, iy)aoo. 

Therefore the matrix C^^V)"^ is also diagonal. To zeroth order. 
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(jxjy I ((2)V) ^ I nx,ny) = (^(jx,nx)(5(jy,ny)/aoo. (45) 

Next, I include first order contributions to *^^V. In the sum over 
(nx,ny), eq. (jS]), I can drop all products of matrix elements which 
are second order or higher. Then the matrix elements in eq. (jH]) 
must be diagonal, with each ket index equal to the corresponding 
bra index (ix = nx, nx = jx, etc.), except for one pair, which I will 
call the off-diagonal pair. For example, suppose the off-diagonal pair 
is ix 7^ jx; then the requirement of at most one off-diagonal matrix 
element implies only two values of nx are possible; either ix = nx 
but jx 7^ nx; or the reverse, ix 7^ nx but jx = nx. Any other choice 
of nx would lead to a product of two first order matrix elements, 
therefore a contribution to the sum eq. (jUj) which is second order 
and negligible. Also, if the off diagonal pair is an x-pair (as in our 
example, where ix 7^ jx) then all y dependence must be diagonal to 
avoid second order contributions: iy = ny = jy. 

For a specific example with ix 7^ jx, consider ix = the ml(M) 
state, jx = the original coherent state. Eq. (jH]) becomes 

= S„,.K,p,.;---|(2V-Mna;,---)(nx,---|(2V |ml(M,);---) 

+ {u^J. I (^V-i I ml(M,))(ml(M,) | (^V | ml(M,)). (46) 

The • • • on the first line indicate suppressed labels iy, ny, jy which 
are all equal and do not change. The intermediate state | nx) = 

1 ml{Lx)) does not contribute; matrix elements 

(ml(L) I (2)V I ml(M)) 

turn out to be second order. 

Since three out of four of the matrix elements in eq. (1461) are 
known to the required order, I can solve for the fourth matrix ele- 
ment. 

(u.,p.;iy I (2V-^ I ml(M.);iy) = -atoM/(aoo)'. (47) 

By systematically considering all ix 7^ jx and iy 7^ jy pairs, I can 
determine all first order elements of the inverse. For example, 

{ix;uy,py I (2)v-i I ix;ml{My)) = -a*l^^/(aoo)^ (48) 
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which will be recognized as a relabeled version of eq. (H7I) . 

Eq. (145|) for the diagonal elements of C-^V)"^ is already correct 
to first order; there are no first order corrections to the diagonal 
elements. Proof: the diagonal elements have ix = jx and iy = jy. 
If also ix = nx and iy = ny, then all matrix elements are diagonal 
and we need a first order correction to a diagonal element. There 
are none (see eq. (H5]) ). If (say) ix 7^ nx, then jx 7^ nx also. This 
contribution to the sum over Ux is second order, and may be dropped 
from the sum. Formula eq. fH5l) for the diagonal elements of (*^^V )~^ 
is correct to first order. □ 

'•^•V can have only an approximate inverse. We can estimate the 
values of the parameters when the series begins to diverge and ^^V no 
longer has an approximate inverse, by examining the first order cor- 
rections to ^^V . All matrix elements of C-^V)"^ have at least one 
factor of aoo in the denominator, where aoo is the ^^V eigenvector 
for the original coherent state. This is as expected: if the volume 
eigenvector is small (if the original state is close to the kernel), then 
the inverse does not exist. To get some measure of closeness to the 
kernel, consider the first order terms. They are down by factors of 
(oqi ora*Q)/aoo- These expressions are small, of order l/\/< L >, 
unless aoo is small, of order a/< L > rather than ^J^T^^JI^. aoo 
is order iy^X^X^"siiie^, where sine denotes the sine of the angle 
between the two angular momenta. When this sine becomes small, 
of order l/< L >, *^^V loses its inverse. 

E The square root of E| 

At the beginning of this sectioii I stated that the square of the 
volume operator is a product of E| times (^^V)^ ; and the resultant 



E^ factor was trivial because the basis states are eigenf unctions 

of E^. This is an oversimplification. The literature contains at 
least two recipes for extracting the square root of volume operators. 
At present it does not seem possible to distinguish between the 
two approaches, using general principles; and I need to state which 
approach I am using here. 

When more than three edges terminate at a given vertex, the 
(volume)^ operator grasps each triplet of edges in turn, generating a 
series of amplitudes, one for each triplet. To parallel a terminology 
from classical optics: should one add first, then take the square 
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root? Or take the square root first, then add? I.e. should one add 
the amphtudes from each triplet, then take the square root of the 
magnitude of the result; or should one take the square root of the 
magnitude of each amplitude first, then add. The literature contains 
advocates for both "add first" [51[7| and "take the square root first" 
[HI [H] choices. For a discussion of the distinct regularization schemes 
leading to each choice see reference [7j. Presumably no final choice 
between the schemes can be made until the two choices have been 
checked in applications. 

In this paper I adopt the "add first" choice, for the following (non- 
rigorous) reason. Consider a four-valent vertex Vn such that two of 
the edges meeting at the vertex are tangent, one ingoing and one 
outgoing. Rotate the gauge at the vertex so that the holonomies 
along these edges are pure Sz, with Sz eigenvalues Mj and Mj 
for outgoing and ingoing vertices, respectively. The vertex has two 
holonomies, 

/•n+l pn 

exp[i / Mf Af dx' ] ; exp [i / Mi Af dx' ] . 

Jn Jn—1 

In the "add first" prescription, the contributions to (volume)^ from 
these two holonomies are added, giving a factor of {Mf + Mi)/2. The 
factor of 1/2 comes from the integrations over half a delta function. 
The volume is proportional to the square root 

v/|Mf + Mi|/2. 

In the "square root first" prescription, the corresponding factor 
would be 

Vl Mi I /2 + Vl Mf I /2. 

In the special case Mf = Mi, it is possible to view these two 
holonomies as a single holonomy passing through the vertex f„. 
When this is grasped by the (volume)^ the contribution is propor- 
tional to Mi with no 1/2, and the volume is proportional to ^y\ Mi \. 
This result for one holonomy equals the Mf = Mi limit of the result 
for two different holonomies, only if the "add first" prescription is 
used. 

The planar case resembles the example just given. Each vertex 
has two holonomies, 

/■n+l pn 

exp[i / MfAfdz];exp[i / MiAfdz]. 



ii-i 
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Given the "add first" choice, the eigenvalues of E^ from the two 
holonomies should be added, yielding the eigenvalue 

(7fi:/2)(Mf + Mi). (49) 

The eigenvalue of the volume operator contains the square root of 
the magnitude of the above expression. 

V The commutator [h,V] 

At this point it might seem we are finished. The basic building 
blocks for the Hamiltonian are the holonomy and the volume op- 
erators; formulas of the previous sections should be adequate to 
evaluate most matrix elements of /i*^^/^^ and *^'^V . 

There is one exception, however. Commutators such as 

-{u,p\ (^V \n)(n\ /i^^/^) I u,p) 

are a small difference between two large terms. The contribution 
from the leading term, | n) =| M,p), cancels out of the difference, 
and it is very hard to evaluate this commutator using only the results 
of previous sections. 

Thiemann and Winkler have shown that, in the classical limit, 
the quantum mechanical commutator equals the classical Poisson 
bracket [1]. 

(u, p I [Oi, O2] I u, p)/ih = {Oi, 02}(u, p) (50) 

On the right, one computes the Poisson bracket, treating the oper- 
ators as classical fields, then evaluates the result at the peak values 
for the coherent state. Since the Poisson bracket is no longer a small 
difference between two large terms in general, it can be evaluated 
using the formulas of previous sections. 

[V,h] commutators occur (for example) in the Euclidean part 
of the Hamiltonian, which contains the following terms (in a field 
theoretic formulation, before discretization on a spin network) 



F' (^1 /e 



(3). 



FfyE^z^'s^zBA/e^'^. (51) 
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In this section A,B = X,Y only. In spin network theory for the 
planar case, these terms become commutators 

2FtyTr{aAh^[h-\^^\]}. (52) 

h is short for h^^^'^\ 

Eq. (1521) is not yet in its final spin network form. Each Fij must be 
replaced by a product of four holonomies encircling the four sides of 
the area ij. This is not the place to discuss the construction of these 
holonomy products; however, one issue does need to be addressed. 
In general a product of four spin 1/2 holonomies will contain all 
the Pauli matrices. For instance, the holonomies replacing F^ oz 
could contain a a matrices, A 7^ Z, or even the 2x2 unit matrix. 
Fortunately, there exist several possible holonomy products which 
reduce to the same Fij in a field theory limit. (For example, proceed 
counterclockwise vs. clockwise around the area.) One can choose 
a linear combination of the possibilities which is pure oz or pure 
(Ta, as needed to reproduce the classical results. In what follows I 
will assume this has been done, so that the only traces required are 
those shown in eq. (1521) . 

Also, one must decide how to choose the holonomies hk in eq. (!52|) . 
In spin network theory the Euclidean Hamiltonian has the following 
form. 

HECxe'J'^hijhk[V,hi;i], 

where hij stands for the product of four holonomies which replaces 
the Fij of eq. flS2]) . When this Hamiltonian acts on a spin network, 
the hk and hij are taken to lie along edges of the network. This 
requirement removes some, but not all of the ambiguity. If the 
vertex is n-valent, there could be as many as n-2 edges which do not 
lie in the ij plane; and any one of these (n-2) edges could qualify as 
support for h^. (In the planar case this ambiguity occurs for ij = xy, 
k = z, since at each vertex there are two z holonomies, one entering 
and one leaving.) I adopt the simplest choice for resolving this 
ambiguity, which is to average over the n-2 possibilities [lOj. When 
combining the n-2 amplitudes contributed by each /i^, I adopt the 
same "add first" prescription as was used for the volume operator. 
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In the planar case it will be possible to evaluate the [h~^, (^V] 
commutator exactly, since h^ is an eigenfunction of *^^V ; see the next 
subsection. One can then evaluate this exact quantum expression 
in the classical, large quantum number limit, and compare the limit 
to the result obtained from eq. (150|) . or equivalently compare the 
limit to the classical field theoretic formula for the commutator. The 
comparison will not support any particular prescription for choosing 
the hk, but does clarify some issues regarding eq. (!50|) . 

A Holonomies along z 

As indicated just above, the hz commutator can be done exactly, 
without using the theorem, eq. (!50|) . The nth vertex has two holonomies 
directed along z, one incoming and one outgoing. Therefore the 
wavefunction at vertex n contains the product 

h{z,n) := exp[i / A^ Szdz] exp[i / A^ Szdz] 

Jn—l Jn 

= exp[i / A^ Midz] exp[z / A^ Mf dz]. (53) 

J n—l Jn 

The total wavefunction is h(z,n)| Ux,Px]Uy,Py). When hz [h'^'^, ^^V] 



acts on this wavefunction, the \J \ ^^-15 | factor in ^^V = W | E| *^^te | 
merely factors out and acts on the coherent state factor. The non- 



trivial part of the commutator is determined by the \j \ E^ | factor 
in the commutator acting on the h~^ factor in the commutator and 
the h(z,n) factor in the state. 



It is straightforward to compute the action of \J \ E^ | on h(z,n) 

and /i~^, since the operator E^ is diagonal in the holonomy repre- 
sentation. We need only compute the action of E|, add (or average) 
amplitudes appropriately, and then take a square root. The holon- 
omy h~^ may be either incoming or outgoing. If it is outgoing, then 

[K\W,]h{z,n) 

i-n+l 

= [exp{-i / Af (±1/2) dz}, E|] h{z, n) 

J n 

= {^K/2)h-'h{z,n)[{Mf + M,) - {Mf + M,t 1/2)]. (54) 
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The ± refers to the [±, ±] diagonal element of the matrix h~^. The 
ingoing holonomy gives an identical contribution. This tells us im- 
mediately that the exact calculation throws no light on how to av- 
erage over the various possible hk- We get the same result, whether 
we use ingoing only, use outgoing only, or take the average. 

I have chosen the prescription which averages incoming and out- 
going contributions, then takes the square root of the magnitude. 
Note however the final two parentheses in eq. flM|) should not be 
combined or averaged; they come from two separate applications of 
the volume operator. 

[K\^f^^]h{z,n) = ,/{^K/2)h;'h{z,n)[^\Mf + Mi\ 



Mf + M.lTvf^l (55) 

where | E| | = ?7E| and r] = ±1 is the phase of the eigenvalue 
Mf + Mi. 

From eq. ( l52l) the result eq. (l55ll should be inserted into the trace 



T{Z,z)h{z,n) := Tr{azK[K\ Jt^E^^]} H 



,z,n 



S±{(±l)v/(7«:/2)U/|M^ + M, 



Mf + Mi\ Tv/'^]}hiz,n) 

= ^{^n/2)[-^\Mf + M,\-7]/2 

+ ^\Mf + Mi\ +7]/2\h{z,n). (56) 

Eq. (l56l) is an exact result. In the limit of large quantum numbers 
I Mf + Mj 1^ 1/2, the radicals can be expanded, yielding 

T(Z, z)-^v/(W2y[V( VI Mf + Mi I) ]. (57) 

This approximate result should be compared to the result from clas- 
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sical field theory (or from tlie tfieorem, eq. (150!) 1 



T{Z,z)h{z,n) -^ Tr{az{l)[l-iAtSzdz,J7]E'z]}Hz,n 



Tr{az{l)[Szdz{6/SE^z)J7iE%]}h{z,n 



= {7]/2)^{^K/{2\Mf + M,\)h{z,n). (58) 

Tliis result agrees with the large quantum number limit, eq. ( 1571) . 
This agreement sheds no light on the procedure for combining 



amplitudes, but does shed light on the theorem, eq. fl50l) . The yrjE^ 
operator is not distributive. When acting on a product of two func- 
tions of holonomies, 



'r^EI [fi(h)f2(h)] ^ f2(h) ^r/EI [fi(h)] + fi(h)^r/E| [f2(h)]. 
This can be checked readily by choosing the fi to be eigenvectors 



of E|. In eq. f l50l) suppose one takes Oi = yrjE^, O2 = /i- The 
theorem seems to state that the action of Oi on /i is independent of 
the state on which the commutator acts (is independent of /2). One 
might be tempted to question this result, given the non-distributive 



character of yrjE^- However, the classical limit of the exact quan- 
tum result agrees with the theorem. 

When the eigenvalues | Mj | of | E| | are large, naively one 



might expect [h^'^, \/\ E| | ] to be order the eigenvalue of -\/| W 



z h 



a/I Mi + Mf I; but instead it is of order 1/(2a/| Mj + Mi |), which 
is the derivative of ^y\ Mi + Mf \. In the next subsection we will 
find that this commutator = derivative structure holds also for the 
commutators of transverse holonomies with the volume. 



B Holonomies along x,y 

Now consider commutators involving /ij, i = x,y. As before, "-^V factors 



into *^^V times \J\ E^ |; but now the W | E| | factors out. When act- 
ing on the h(z,n) factor in the wavefunction, it gives a factor of 
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^1 Mf + Mj I, as at eq. (l49ll . We are left with the following trace, 
which acts on the coherent ket. 

T(A, i) = Trace(aA h; [hr^, (^V] ). (59) 

Since the case i = y may be obtained from results for i = x by a 
relabeling, it suffices to compute this trace for i = x. 

The theorem, eq. fl50l) . gives for the commutator in T(A,i) 

= (7«:/2)[5B,/.;i)] + [CE^ce^''/2('V], (60) 



where ^^V is y C E^ -^c- C'^'" and ^ = ±1 is the sign of the eigenvalue 

of ((^V)^ . The anticommutator arises on the second line because 
the holonomy is a loop, beginning and ending at the vertex where 
the E acts. Therefore the grasp occurs on the far left and right of 
the holonomy. 

The *^^ V in the denominator of the formula eq. fl60|) requires some 
discussion. Even though *^^Vhas vanishing eigenvalues, there is no 
divergence. The "^^V operator will be acting on coherent kets which 
are approximate eigenvectors of ^^Vwith large eigenvalue. For this 
situation, from the discussion of (*^^V)~^ in section [TV] the zero 
eigenvalues may be neglected. '^^V connects the original large eigen- 
value state primarily to other large eigenvalue states, rather than 
states in the kernel. 

The ^^V factor also does not give rise to a factor ordering ambigu- 
ity, even though it does not commute with the h~^. The difference 
between two orderings equals a commutator, which is small. 

h,{l/m) = (l/(2V)/i. + [/^,,l/(2V]; 

= order (/i,/(2V)([/i,,(2)v]/(2V). (61) 

On the second line, I assume h^ is outgoing, so that the *^^ V operator 
overlaps with h^ on the left side of hx- An ingoing holonomy would 
give the same final order of magnitude. In the present case, as in the 
previous subsection, the commutator resembles a derivative (now a 
derivative with respect to L rather than M) in that the commutator 
lowers the power of L by one. To see this, note ^^V is order -yjE^Ly, 
from the two E operators in *^^V . From eq. ( l60l) . [hx, "-^V] is order 
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Lyl ^L^Ly, down by a derivative with respect to L^- (The Sb is just 
a Pauh matrix over 2; it is not order L.) Therefore the two order ings 
in eq. (|6T1) differ by a factor which is suppressed by [h^, (^V]/*^^V = 
order l/Lx- In the classical limit, factor ordering is not a problem. 
1 now substitute eq. fl60l) into the expression for the trace, eq. fl59l) . 

= [5ab + /^('^ {hx)AB] (7/^/2) C esc Ey (2V ; 
D«(/i) = DW(-0 + 7r/2,^,0-7r/2). (62) 

The second line uses the rotation property of the su(2) generators, 
eq. (jHSD- The matrix on the last line is a full-angle rotation, not a 
half-angle rotation; we are on the SU(2) side of the SU(2) vs. 0(3) 
divide discussed in section II of paper 1. 

Eq. (162!) is correct for T(A,x). For T(A,y) replace hx^hy and 

A note on commutator = derivative: taking the commutator does 
reduce the magnitude in the same way as taking a derivative; but 
the commutator does not preserve directions. Compare matrix el- 
ements of (^Vto those of [h-\ (^V], eq. ([60]). I replace the SU(2) 
objects /i^^/^) by 0(3) objects D(1)q^ to facilitate order of magnitude 
comparisons: 



< ('V > ~ Y^< EJ, X E^ > eMN] 

<[^(1U^'V]> ~ < D(l)o^ X E^ > e,W < (^V > . 

The net effect is to replace < Efj > by < D(l)g^ >. From eq. (I9l), 
< E^ > is order < L^ > 

< EIj >=< Lx > (hncos/i -n X V)Msinyu), 

while its replacement is order unity. Evidently the commutator 
reduces magnitude in the same manner as a derivative. 

However, < EJ^ > also contains a unit vector, and information 
about this vector has been lost. Presumably the commutator = 

derivative relation was exact for commutators with -\/E|:, because 
only one direction is involved. For other operators, the relation is 
useful when one is counting powers of L. 
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VI Conclusion 

One technique used the present, planar calculation depends only 
on angular momentum theory, therefore applies also to the general 
SU(2) case (or to any symmetry which has recoupling coefficients 
analogous to 3 J symbols). Given a coherent state with holonomy 
peaked at a matrix u, all peak values should be functions of u. Since 
the basic planar operators are spin 1, some step in the calculation 
should produce low-spin representations D^^^ of u. In this paper I 
have used the formulas for the rotation behavior of Clebsch-Gordan 
coefficients, appendix \^ to reveal this u dependence . The basic 
operators in the general SU(2) case are spin 1/2; therefore the same 
formulas will produce low-spin representations D^^/^) of u. 

Those formulas also lead to a clean separation between leading 
terms and the small correction terms. Study of the SC terms reveals 
that they emphasize small fluctuations, values of the basic variables 
M and L which are near, but not at the mean values < M > and 
< L >. Further, the standard deviation parameter t cannot be 
taken too different from l/< L >, or these SC terms will become 
large. 

The coherent state approach is not like perturbation theory. One 
cannot improve on the initial approximation by taking more terms, 
essentially because perturbation theory uses a complete set, whereas 
coherent states are overcomplete. 

Can one construct a complete subset of the overcomplete set, a 
subset suitable for perturbation theory? In effect, I did this in sec- 
tion IIVI The subset consists of the original coherent state, which 
has Gaussian distributions in L and M, plus SC states with distri- 
butions given by moments of the Gaussians. The subset was used to 
construct a perturbation expansion for the inverse of *^^V . Indeed, 
all the formulas giving the action of the holonomy and E operators 
on the coherent state may be considered as perturbation expansions 
in this complete subset. 

The perturbation calculation in section IIVI certainly is helpful. 
Using the perturbation expansion, one can put quantitative limits 
on when *^^V no longer possesses an approximate inverse. 

Also, the technique works quite generally. One can use moments 
to construct a complete subset, for calculations involving any kind 
of coherent state. 
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Is this perturbation approach useful for extending coherent state 
calculations to smaller values of L? Since higher terms in the per- 
turbation series fall off as 1/VX, perhaps one could extend the cal- 
culation to low values of L, L < 100, by including higher moment 
terms in the series. 

However, there is a possible problem with spreading of the coher- 
ent wave packet. Coherent states are useful because their eigenval- 
ues are strongly peaked at one central value. If this peak broadens 
rapidly in time, the states lose their coherent character. In a follow- 
on paper, I estimate the rate of spreading of a coherent state, and 
for small L the rate of spreading appears to be quite rapid [2] . The 
perturbation expansion is certainly useful in many contexts, but it 
may not be well suited for reaching small values of L. 

Future work should choose a specific Hamiltonian, then check for 
wave propagation and correct commutation relations in the coherent 
state limit. Although I have not written down a specific Hamiltonian 
in this series of papers, I have listed constraints which limit the form 
of the Hamiltonian. (See for example the remarks following eq. (!52l) .) 

A Identities Involving D and S 

For the convenience of the reader, this appendix includes brief deriva- 
tions of two well-known identities involving the rotation matrices D 
and generators S. Throughout, I do not use the relation D~^ = D], 
so that the formulas in this appendix remain valid for matrices in 
SL(2,C). My conventions for Clebsch-Gordan coefficients and rota- 
tion matrices are those of Edmonds [TT] 

The first identity relates matrix elements of S to a Clebsch- 
Gordan coefficient. 



(L, M' I Sa I L, M) = VL(L + 1)(L, M' | L, M; 1, A), (63) 

where the Clebsch-Gordan coefficients for LI (^ L2 = L3 are writ- 
ten in a braket notation as (L3, M3 | LI, Ml; L2, M2). The M 
dependence of the right hand side of eq. ( 163|) is required by the 
Wigner-Eckart theorem applied to a spin one operator. To check 
the scalar coefficient \/L{Ij + 1), square both sides and sum over M 
and A. 
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(L,M' I S^ I L,M') = L(L + 1)(L,M' | L,M'). 

The phase of the scalar coefficient can be verified by checking a 
simple example. 

The second identity is the formula for reducing a product of two 
rotation matrices. (I suppress the labels Lj, which are obvious from 
context.) 

DmMiDN2M2 = J2{N3\ N1N2)Dn3M3{M3 \ M1M2); 

L3 

{N3 I mN2)DmMi = Dm3M3{M3 \ MlM2)D];l^j,^. (64) 

The second line is a rewritten version of the first. 

As an illustration of these rules, I obtain the usual rotation prop- 
erty for the Sa vector operator. Insert eq. (|63|) for S into the second 
line of eq. (I64l) . and restore the L's. 

(L, N3\Sa\ L, Nl)D^j}l,, = Df,^j,{L, M3 \ Sb \ L, Ml)D^i\-' 

SaD(^) = D^^^Sfil^gl"'. (65) 

The last line uses a matrix notation to hide some indices. This 
formula is also valid for matrices D in SL(2,C). 

For components of vectors, this paper use both Cartesian (X,Y,Z 
or 1,2,3) and spherical (+1,-1,0) indices. Strictly speaking, the two 
should be treated slightly differently when forming dot products. 



A-B = AxBx + AyBy + AzBz 

= (A_)*fi_ + (A+)*fi+ + A353. (66) 

Since spherical components are complex, the second dot product 
resembles the dot product in SU(3) rather than 0(3). For the most 
part I have omitted the complex conjugation stars, trusting to the 
reader to know when to put them in. Note in most (but not all) 
cases, no star is needed because one of the indices refers to a final 
state, which can be considered starred. For example, there should 
be a star on one of the factors in eq. flTSl) : but no star is needed in 
the sums over magnetic quantum numbers in eq. fITU]) . 
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B The matrix D^^\u) 

This appendix calculates the rows of D^^\u), which form a natural 
basis for the vector operators E and h. More precisely, from eq. (1131) . 
the relevant basis vectors are 

exp[i /? B] D W(u)ba; A, B = +1,0,-1. 

The subscripts (B,A) give the components A of unit vector B. 
One may evaluate the components of each vector, starting from 



D(i)(u)ba = exp[i(/? - 7r/2)(A - B)]S'\a/2) 



BA, 



with 



dW(«/2) 



(l + cos)/2 +sin/v^ 
— sin /^/2 cos 

(l-cos)/2 -sin/v^ 



(1 - cos)/2 

+ sin / V2 

(1 + cos)/2 



(67) 



cos = cos(a/2),sin = sin(a/2); rows and columns are labeled with 
spherical components in the order (+1, 0, -1). 

From eq. ([8]), the basis vector B = is just the unit vector V^(u), 
i. e. V^(h) with h = u. This vector has spherical and Cartesian 
components 



V{u)a = D^'\u)oa 

= {t sin(a/2) exp[±i{(3 - 7t/2)]/V2, cos(a/2)) 

= (sin(a/2) sin(/5), - sin(a/2) cos(/?), cos(a/2)). (68) 

Spherical components are listed in the order (±,0); Cartesian com- 
ponents in the order (X,Y,Z). 

As for the rows i? = ±1, it is best to form linear combinations in 
order to get real vectors. One linear combination turns out to be n, 
the axis of rotation for u, while the orthogonal linear combination 
turns out to be n x V, the third axis of an orthogonal coordinate 
system (n, V,h x V). 
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V2nA = -exp[+ip]D^'\u)+i,A + exp[~ip]D^'\u).i,A 

= -^(rfS,A + rf-M)expM(/5-7r/2)] 

= {Texp[±t(3],0) 

= y2(cos(/?),sin(/?),0); (69) 

i\/2(n X V)a = -exp[+tf3]D^^\u)+i,A - exp[-i(3]D^^\u)-i,A 

= -i{d^liA - ^-m) exp[tA{f3 - 7r/2)] 

= i{Tcos{a/2) exp[±i{f3 - 7t/2)], -V2sm{a/2)) 

= iV2{cos{a/2) cos{f3 - 7r/2), cos(a/2) sin(/? - 7r/2), - sm(«/2)). 

(70) 

For convenience I also record the inverses of eqs. (169!) and (170|) . 

l^^'Hw)oA = V{u)a; 
exp[+i(3]D^^\u)+i,A = {-fiA - i{n x V)a)/V2; 
exp[-z/3]D«(«)-i,A = {hA-tinxV)A)/V2. (71) 

The vector V^(u) will turn out to be the peak value of V^(h); 
therefore (since the relation between V and h^^^'^^ is invertible) u^^^"^^ 
will be the peak value of h^^'"^' . We could have used the notation f 
for the unit vector which we labeled V^(u), because the peak value of 
radius (for the fictional electron moving in a central force) is given 
by the vector V^(u). 

C The matrix L'(^)(H) 

This appendix derives the properties of the Hermitean factors D(H) 
occurring in the coherent state. The initial formulas will be valid for 
general p and p; later results will use the assumptions p = (cos(/3 + 
/i), sin(/9 + /i), 0); i. e. pa = 0; and e^ ^ 1. 

Most of the mathematical techniques used in this appendix are 
a direct steal from Thiemann- Winkler paper II [12], with one sig- 
nificant exception: I make no use of traces. Because Thiemann 
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and Winkler deal with the general SU(2) expansion (matrices Dmn, 
both M and N summed over) they are able to recast their results for 
D(H) as theorems about class invariants, the traces Dmm ■ In niy 
case the expansion matrices are Dqm, sum over M only; I have not 
been able to recast my results as theorems about traces. Instead, 
in order to obtain manageable forms for D(H), I use the assump- 
tion e^ ^ 1. This assumption is not a serious limitation unless one 
wishes to extend the calculation to very low values of L of order 10. 
(For discussion of this point, see the estimates given for the size of 
t, in appendix [D1 These estimates in effect also limit the size of p.) 
To obtain an approximate form for D(H), I write down the power 
series expansion for D(H), then use the e^ 3> 1 assumption to sim- 
plify the series. Start from D{H)'^^^'^\ which has the form 



D{H) 



1/2 



exp[p- S] 

cosh(p/2) + p3 sinh(p/2) sinh(|3/2)(pi — ip2) 
sinh(p/2)(pi + ip2) cosh(p/2) — p^ sinh(p/2) 

b 



a 
c 



d 



(72) 



(Despite the superscript 1/2, this is a matrix in 0(3).) Given the 
abed, the D(H)^^^ for arbitrary L may be computed from the series 



D^'^\H)mn = 

{adf{a/h)^{h/df^^{L + M)\{L - M)\{L + N)\{L - N)\ 

{bc/adY 



'{L-M -k)\{L + N -k)\{M -N + k)\k\ ' 



(73) 



This formula is valid for H an arbitrary complex matrix in SL(2,C), 
not just for H a Hermitean matrix. 

A Large p limit of D(H) 

We now assume e^ large; pa may be large, but not too close to ±1: 
(1- I P3 I) > 2e'P.) Then 
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cosh(p/2) = sinh(p/2) 

^ 6^/72; 
ad = be; 

D^'^\H)mn = (ad)^(a/6)^(6/rf)^V(^ + M)\{L - M)\{L + N)\{L - A^)! 
■Efep - M - k)\{L + N- k)\{M -N + k)W]-^ . (74) 

The series in k can be summed using the addition theorem for bi- 
nomial coefficients: 

The eqs. (EH) and ([75]) give 

D^'^\H)mn = {adf{a/b)''{b/d)^ ^ ^^^^' 

^ ' ^ ; V / ; V / ; ^(L + M)\{L - M)\{L + N)\{L - N)\ 

= (Pi - zp2)'''-'^(exp(p/2)/2)2^[l +p3]^+^[l -pa]^-"" 

(2L)! 



v/(L + M)\{L - M)\{L + N)\{L - N)\ ' 

At this point one can prove: let M and N denote the peak values 
of M and N, i.e. the values which maximize | D{H) \. Then 

M/L = N/L ^ p3. (77) 

Proof: to find the peak value of (say) N, compute the first difference 
of the square magnitude of the N dependence of D(H), and set this 
ffist difference equal to zero. 

^^'V(iV) := f{N+l)-i\N); 

f{N) = \a/b\''' /[{L + N)\{L-N)\]. 
= ,5«/(Ar) 

, /, ,2 L-N 

oc a/b = 1: 

' ' ' L + A^ + l 

I a P - I 6 |2 
I a I + I I 
N/L ^ p3. (78) 
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(76) 



On the last line I have used the values of abed from eq. (1721) . The 
proof for M is identical except for the replacements (a,b) ^(b,d). 

D 

B Small p3 limit of D(H) 

In the main body of the text I focus on the case pa = 0. From 
eq. fl77|) of the last subsection, in this limit the important values 
of M and N satsify L » M,N. Therefore one can use Stirling's 
approximation for the factorials in eq. fl76l) . for example 



(2L)! (2L2L 



(L + M)!(L - M)! 0F(L - M)l-m+i/2(l + m)l+m+i/2- 
Now use 

(l + x/ri)"" = exp[nln{l + x/n)] 

= exp[x-xV2n + ...]. (80) 

Take n = L + 1/2, x = ±M. Also, write L^^ as (L + 1/2 - 1/2)^^ 
and apply eq. ( IHOl) to this factor. 

(2L)! ^ (2^i-exp[-MV(L + l/2)] 

(L + M)!(L-M)! v/7r(L + 1/2) ^ ' 

To obtain a result valid near p^ = 0, as well as at p^ = 0, assume 



Ps < order 1/a/L + 1/2. Then apply eq. flHOj) to the [1 ±^3] factors, 
with now n = L+1/2, x = ±p3(L + 1/2) . For example, 

[1 + p,]L+^ = [1+ p,{L + 1/2) /{L + l/2)](^+l/2)[l+(iV-l/2)/(L+l/2)] 

^ exp[p3(^ + l/2)+P3(iV-l/2) 

-(p3)^(^ + l/2)/2], (82) 

and similarly for the [1 — ^3] factor. Inserting eqs. ( IHTl) and (1821) into 
eq. (1761) yields 

^iL) ~ /^ „-^ \M-N exp(pL) 



D{HY^>^ = {p,-zp. 



^Ti{L + l/2) 
■ exp{-[M - p,{L + 1/2)] V2(L + 1/2)} 
• exp{-[iV - ps{L + 1/2)] V2(L + 1/2)} (83) 
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The M (and N) dependence of D(H) is peaked at M = p^{L + l/2), 
with the squared width of the Gaussian equal to \/L{L + 1) = (L + 
1/2). This is aheady a bit more than we need for the main body of 
the paper. 

Eq. fl55]) demonstrates Gaussian behavior in N and M. To obtain 
Gaussian behavior in L, multiply D(H) by the other exponential 
factor in the coherent state. 



exp(-tL(L + l)/2)D{H)% ^ exp(-tL(L + l)/2)e^P ■ ■ ■ 
= exp[-t((L + 1/2) - j9/t)V2 + py{2t) + t/8 - p/2] ■■■ . (84) 

The ■ ■ ■ indicates irrelevant factors which are bounded for large L. 
On the last line one can neglect exp(t/8) = 1. Eq. fl5^ is a Gaussian 
in L with mean < L + 1/2 >= p/t and standard deviation 1/Vt. 

Since t is small, the standard deviation is very large. However, 
what counts is ai/ < L + 1/2 >= \/t/p, which is small as required. 

When p3 = 0, at eq. (jTj) I have introduced the notation 

p = (cos(/9 + /i), sin(/5 + /i), 0) 
for p. For this value of p, D(H) becomes 

exp(-tL(L + l)/2)D{HY^i - exp[-t((L + 1/2) - p/tf/2] 
■ exp[-MV2(L + 1/2)] exp[-Ny2{L + 1/2)] [l/y/n{L + 1/2)] 

■exp[pV(2t) -p/2](expH(/5 + /i)])^^-^. (85) 

D Small Correction (SC) Terms 

This appendix discusses the nature of the small corrections SC, 

operator | coh state) = <operator> | coh state) + SC. 

The first subsection constructs a set of states which are orthonormal, 
and uses them to expand the SC for the E operator. The coefficients 
multiplying these states are shown to be suppressed by factors in- 
volving the small parameters l/\/< L > and \/i. A second section 
shows that the SC terms emphasize values of the parameter near. 
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but not at the average value. The term "near coherent" is intro- 
duced to describe this behavior. A third subsection introduces a set 
of states appropriate for expanding the SC terms for the holonomy 
operator. A final subsection argues that the parameter t should be 
taken to be order 1/ < L > in order to minimize the size of the SC 
terms. If t is replaced by a number of order 1/ < L >, everywhere 
in the factors multiplying the SC terms, then all the SC terms turn 
out to be suppressed by factors of the same order, 1/V< L >. 

A SC States for the E Operator 

At eq. (ITBIl I replaced 

L^ < L > +AL; M^AM, 

then asserted that the terms proportional to < L > represented the 
dominant contribution. I must now examine the terms involving 
AX, X = L or M, and show that they are small. 

From the previous appendix, eq. fl85|) . the original coherent state 
is proportional to Gaussian factors coming from the D(H) factor. 
Therefore the SC terms are proportional to the first moments of 
these Gaussians. 



u,p) ^ Y,D^''\hu]\ 



)0M 
L,M 



■ exp[-t((L + 1/2) - p/t)V2] exp[-MV2(L + 1/2)]; 
SC terms oc ^D(^)(/iMt)oAf [ALor AM] 

L,M 

■ exp[-t((L + 1/2) - p/t)V2] exp[-MV2(L + 1/2)]. 

It is difficult to carry out the sums over L and M, because the 
D^^\hu\)oM factor is difficult to approximate and may have com- 
plicated L and M dependence. 

There is a simpler way to estimate the order of magnitude of the 
SC terms, without knowing in detail the M and L dependence of 
D^^\hu\)oM- For the SC terms involving M and AL, define the 
states 

I ml(M)) := Ar(ml(M))^((2L + l)/47r)exp[-tL(L + l)/2] 

L,M 

■[D^'^\hu])oMMD^^\H)]Mo; (86) 
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ml(L)) := Ar(ml(L))^[(2L + l)/47r]exp[-tL(L + l)/2] 

L,M 

■[D^''\hu])oMALD^'^\H)]Mo. (87) 



The notation mp(X) denotes the pth moment of the variable X. The 
above states are identical to the original coherent state | u, p) except 
for a different normalization factor, 

N^N(ml(X)), 

and one power of AX in the summand. In terms of these states, 
eq. (fT5|) becomes 



(7K/2) -^E^ \ u,p) = < L > {fiA cos fi — h X V) A sin fi) \ u,p) 

+ {N/N{ml{L))) I ml(L)) 

+ {N/N{ml{M))[VA + i sin fifiA 

+i cos fi{nxV) a] \ml{M)). (88) 

Evidently this replaces the problem of evaluating the sums over L,M 
by the problem of determining the normalization ratios N/N (ml (X)). 
This may seem like replacing Tweedledum by Tweedledee, except 
the dangerous factors of D^^\hu])QM drop out when calculating 
norms. 

I now prove the following: 

1 = N^ exp[p^/t -p]^y<L + 1/2 >/2TrVt; 
N/N{ml{M)) ^ V< L + 1/2 > /2 : 

N/N{ml{L)) ^ v^VM- (89) 

I begin with the first line of eq. fl89l) . Orthogonality for the D(h) 



IS 



f D(S(h)D2(h)* = 5l,l'5m,m'(2L + l)/47r. 
Jn(h) 



'C(h) 

From this and eqs. ([1]) and 
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1 = {u,p\u,p) 

= N'Y,^xp[-tL{L + 1)][(2L + l)/47r]DS(<7t)*Z}lS(^t) 

L,M 

= N^ Y^ exp[-tL(L + l)][(2L + l)/47r] 

L,M,N,N' 

= A^2^[(2L + l)/47r]exp[-tL(L + l)] 

L,M 

■D^iiAH)DfkH). (90) 

The u and h dependence have disappeared. 

The next step is to carry out the sum over M. Compare eqs. (jl]) 
and ([3]): when going from h^/'^ to H'^^^ we make the replacements 

im ■ S^^S ■ p. 

I.e. replace magnitude and direction as follows. 

i 6'-^p; (p-^P + yU (-v^ hi— >p). 

For the angles, see eqs. ([5]) and ([7]). Therefore the Euler decomposi- 
tion of D(H) follows from the Euler decomposition of D(h), eq. (HI). 

D(h)(-0 + 7r/2,^,0-7r/2)^D(H)(-/?-/i + 7r/2,-ip,/? + /i-7r/2). 

Therefore the D's on the last line of eq. (I90ll equal 

D(L)(-/5 - /i + 7r/2, -2ip, /3 + /i - 7r/2)oo. 

Approximate this factor using eq. (|85|) with p— ^2p, t^2t. 

1 ^ N^ V[(2L + .)..,]exp[-t((L + l/2)-p/t)^ + pV(t)-p] 
^ vML + 1/2) 

(91) 
Replace the sum over L by an integral: 



^(AL = 1) = (1/v^)5^A(v^(L + 1/2)-p/v^:=m) 

L L 

= (1/Vt) f du. (92) 
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Elsewhere in the integral, replace L + 1/2 by its peak value 

<L + l/2 >=p/t. 



Eq. (1911) then gives the first line of eq. 

The remaining two lines of eq. fl89l) may be proved using similar 
approximations, with one exception. The calculation of N(ml(M)) 
resembles the calculation of N, eq. fl9T|) . except for an additional 
factor of M^ in the summand, so that the sum over M cannot be 
carried out immediately. Instead, the sum over M may be replaced 
by an integral over a variable w, using 

Sm(AM = 1) = v^L + 1/2 J:{AM/^/L + 1/2 := Aw) 

^ ^L + 1/2 / dw. (93) 

D 

Eq. (I89ll implies that the SC terms are suppressed: from eq. ( l88l) . 
the leading term is order < L >; therefore the SC terms are down 
by factors of order 



N/N{ml{M)) <L> = 1/V2 < L >; 
N/N{ml{L)) <L> = l/{V2i<L>). 

At first glance the formula for N/N(ml(L)) looks dangerous, 
because of the small factor of \/i in the denominator. However, 
t = p/ < L + 1/2 >, from eq. (jHSD- Therefore 

N/N(ml(L))/ < L >= l/v/2p < L >. 

Since I am taking e^ ^ 1, both SC terms are down by at least 
1/V< L >. 

Section [IV] constructs a series expansion for the inverse of the 
volume operator and determines the values of the parameters where 
the series diverges. To evaluate the terms in the series I need the fol- 
lowing theorem, which describes the action of E on the first moment 
states. 

E^ I ml(M)) = < E^ >| ml(M)) 

+b I m2(M)) +c\u,p} 
+d\ml{M),ml{L)), (94) 
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plus a similar theorem for M ^h> L, where 

b,d = order [N{ml{M)/N{b,d)]; 

c = order [N{ml{M)/N] <L> . (95) 

N(b,d) denotes the norm of the state with coefficient b, d. 

In words, eqs. ( !M|) and (p5ll imply that | ml(M)) is an approxi- 
mate eigenfunction of E , with the same eigenvalue as | u,p). Also, 
the normalization ratios in b, c, d can be evaluated using the same 
sum-becomes-integral techniques that establish eq. (15^ . The quan- 
tities b, c, d turn out to be order V<L >. Therefore the b, c, d 
terms are small corrections, since < E^ > is order < L > from 
eq. 0. 

The proof of eqs. flMl) and fl95|) for the most part uses the same 
techniques as the proof for the action of E on the original ket, 
eqs. (|88|) and fl89|) . and there is no point in giving the proof in detail. 
However, a qualitative sketch of the steps in the proof may be in 
order, since the c term in eq. fl9^ is perhaps unexpected. When E 
acts on the original ket, it gives back the original ket, or adds a fac- 
tor of M. In moment notation, the original ket is | r7iO(M), mO{L)). 
The presence of a c term implies E can also remove one power of M. 

I sketch the proof of eqs. (19^ and (p5!) . It is a repeat of the proof 
for the action of E on the original ket, eqs. (IHH!) and (!89|) . except for 
the following changes. In section [TTl eq. (TTOl) . relabel on the left: 

I u,p)^ I ml(M)) (96) 

On the right, replace the factor of D(H) in the ket by 5*^0(11) (a 
concise way to replace D{H)mo by M D{H)mo)- At the next step in 
the derivation, eq. (ITT]) , the M becomes M - B, due to a relabeling M 
— >^ M - B. Therefore in eq. (ITl]) . (relabel as at eq. (1^ . and) replace 
the C's. 

C(L,B)^C(L,B)(M-B). 

The derivation now moves to this appendix, eq. flHHj) . Again, 
relabel on the left, using eq. fl96|) . In the kets on the right, insert 
an extra (M - B) into the summand for each ket. The M term (in 
M - B) changes every ket on the right in eq. flHHj) to a ket with one 
higher moment. 
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1 u,p) - 


> 1 ml(M)); 


nl{M)) - 


> |m2(M)); 


ml(L)) - 


> |ml(M),ml(L)) 



(This change is accompanied by a corresponding relabehng of the 
norms N(i)). The M term therefore yields the leading term in 
eq. ( 19^ . plus the b and d terms. 

Next, consider the effect of the -B term in (M - B). Return to 
eq. flT^ for a moment. We are now inserting factors of C(L,M)(-B) 
into the summand. The -B is just =i=l or 0, independent of M, while 
the expansion of C(L,M) starts off with a leading < L > term (plus 
smaller corrections proportional to AX, X = L or M). This leading 
term no longer contains any factor of M. We have removed a factor 
of M, and returned to the original coherent state. 

The factor of < L > (in the expansion of C(L,B)) explains the 
factor of < L > in the order of magnitude given for the coefficient c, 
eq. 0951) . The remaining factor in c, the ratio of norms, is given by 
eq. (159]) . Despite the factor of < L >, c is only order V< L >, which 
makes it a small correction to the leading term, of order < L >. The 
remaining ratios of norms in eq. fl95l) may be worked out in the same 
manner as the ratios for the original coherent state. □ 

B Near Coherent States 

The states | ml(X)) emphasize values of X which are near, but 
not at the average value < X >. To see this, note the new state 
I ml(M)) (for example) differs from the original coherent state by 
the replacement 

~ Mexp[-MV(2(T^)]. 

D(H) is Gaussian in M with a peak at M = 0; see eq. (1851) . Therefore 
the above function is an odd function of M with a zero (rather than 
a maximum) at M = 0, a peak located at M = +o"m, and a valley 
located at M = -aM- (From eq. ([HS]), aM = a/L + 1/2.) 

The shape of this curve suggests that the new coherent state is 
the difference of two simpler coherent states. In fact. 
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Mexp[-MV(2aI,)] ^ (aMv^/2) {exp[-(M - aM)V(2aM)] 

_exp[-(M + aM)V(2alf)]}- (97) 

(To prove eq. (!97Il . multiply out the exponents on the left, and ex- 
tract common factors. The expression will contain an overall factor 
of sinh(M/aM)- Expand this, neglecting terms of order 

(M/aM)=^ = (M/v/L + 1/2)3. 

These terms are small, because the Gaussians suppress terms with 

{M/aM) > 1. □ 

The state | ml(M)) has an M distribution which is the difference 
between two very similar Gaussian distributions. This explains why 
its contribution is so small. 

In the main body of this paper I focus on the case p^ = 0. How- 
ever, in appendix O I discuss briefly coherent states which have 
small ps, Lp3 < ordervT . When p^ is non-zero, the Gaussian in M 
is peaked at < M >= p3{L+l/2) rather than < M >= 0. Therefore 
the Gaussians in eq. (!97|) correspond to small, non-zero values of ps: 



P3 = ±1/a/L + 1/2. The M distribution of | ml(M)) emphasizes 
small fluctuations in ps, near but not at the original mean value 

P3 = 0. 

Similar remarks could be made about the state | ml{L)). It is the 
difference between two coherent states having peak < L > shifted 
from < L > to < L > ±aL =< L > ±l/^/i. Since < L + 1/2 >= 
p/t, these states emphasize small fluctuations in the magnitude p, 
rather than fluctuations in a component, p^. 

I will refer to states such as | ml{M)) and | ml{L)) as "near- 
coherent" states. Whereas the distributions of the coherent state 
establish a peak value for each variables, the distributions of near- 
coherent states emphasize the fluctuations near the peak value. 

In a certain sense the near-coherent states are coherent. If a 
holonomy or E operator is applied to the near-coherent state, the 
near coherent state is found to be an approximate eigenvector with 
the same approximate eigenvalue as the coherent state. Even though 
the distribution for the near-coherent state may have a zero where 
the coherent state has a peak, the eigenvalue (= average value) is un- 
changed, because the off-center Gaussian peaks of the near-coherent 



48 



state are symmetrically distributed around the original peak value. 
This is the qualitative explanation for the < E^ > term in eq. (IMI) . 

C SC States for the Holonomy 

As with the E operators, I estimate the order of magnitude of the 
SC terms by calculating appropriate norms. Since the holonomy 
produces states containing D^^\h) with L' = L ± 1, the states 
I ■ml{X)) are not enough, and I will need the following additional 

states: 

|L±) := iV(L±)^((2L + l)/47r)exp[-tL(L + l)/2] 

L,M 

|L + l,ml(M)) := A^(L + l,ml(M)) ^((2L + l)/47r) exp[-tL(L + l)/2] 

■D^''^'\hu])oMMD^''\H)Mo- (98) 

In terms of these states, eq. (128!) becomes 

D«(/i)o,A I u,p) ^ N/N{L+) I L^)Va 

+N/N{L + 1, ml(M)) I L + 1, ml(M)) 
•[— cos fifiA + sin/i (n x V")^]/ < 2L + 1) > 
—iN/N{L_) I L_)[sin/inA + cos/i(r?, x V)a\ 
+order(M/L)^l/L). (99) 

The order of magnitude of the SC terms may be determined by 
computing the following norms: 

N/N{L+) ^ 1; 
N/N{L^) = ^/t/2; 



N/N{L + l,ml{M)) ^ v^<L + l/2>/2. (100) 

Proof: the last line of eq. fllOOl) is easiest to establish. Because 
I L+l,ml(M)) differs from | r7il(M)) only in the replacement of 
D^^^^\hu\) by D^^\hu\), and the D{hu\) factors drop out anyway 
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when computing norms, N/N(L+l,ml(M)) is the same as N/N(ml(M)), 
eq. dHH]). 

To determine the N(L-i-), rewrite the L± states as follows. 

|L±) = iV(L±)5^((2L + l)/47r)exp[-tL(L + l)/2] 

L,M 
L,M 

•(1/2) ^[(±1) exp{-t[L T 1 + 1/2 - p/tf/2}] 
± 

exp{-MV[2(L + l/2)]} 



• exp{p72t - p/2 - iM{(3 + /i)}- 
iV(L±) ^((2L + l)/47r)Z}(^)(/int) 



^T^{L + 1/2) 



OM 



L,J\/ 



•exp{-t[L + l/2-p/t]V2} 
■ exp{pV2t - p/2 - iM{l3 + /i)} 



cosh[t(L + 1/2 - p/t)] 
sinh[t(L + 1/2 - p/t)] 

exp{-MV[2(L + l/2)]} 



^AL + 1/2) 



iV(L±)^((2L + l)/47r)D(^)(M) 



OM 



L.Af 



•exp[-tL(L + l)/2]Z}(^)(i7)Mo 
cosh[t(L + l/2-p/t)] 
sinh[t(L + l/2-p/t)] 



(101) 



On the third line I have relabeled L ± 1 = L, used eq. ( 1851) to re- 
place the D(H) by Gaussians, and then dropped the tildes. From 
the last line, the | L±) states are just the original states times an 
additional factor of cosh or sinh. When the state is squared to 
determine a norm, this factor becomes cosh^[v^u] or sinh^[-\/tn] 
(As at eq. fl92|) . the variable L is replaced by a variable u, and 
the sum over L is replaced by an integral over u .) For the L_|_ 
state, cosh^[-\/tM] = 1. The normalization integral for N(L+) re- 
duces to the normalization integral for N, and we get the first line 
of eq. (IIOOI) . For the L„ state, sinh^[v^M] = tv?. This normaliza- 
tion integral should be compared to the normalization integral for 
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I ml{L)). That integral contains a (AL)^ = u'^/t factor. Therefore 
in {N/N{ml{L)))'^ = l/(2t), eq. (l89l) . move the t from denominator 
to numerator to get {N/N{L_)f = t/2. D 

The states | L±) introduced at eq. (^5]) are superpositions of the 
states I L ± 1) studied in section [XTTl 

|L±) = (|L + l)±|L-l))/2. 

In section UTT] it was shown that the states | L it 1) are identical to 
the original coherent state, except that peak angular momentum is 
shifted by one unit, from < L > {= p/t-1/2) to < L > ±1. This 
shift by one unit is very small compared to the width of the Gaussian 
peak (1 -C (Jl = l/y/t). Therefore the superposition which adds 
the two Gaussians, | Iv+), closely resembles the original coherent 
state. The superposition involving the difference, | L^), on the other 
hand, emphasizes fluctuations near the peak in L, rather than the 
peak. This state is near-coherent, very similar to the state | ml{L)). 
Consistent with this classification into coherent vs. near-coherent, 
the normalizations for | L±), eq. fllOOp . make | L+) dominant and 
I L_) a SC. 

D Estimates of the parameter t 

It is difficult to put significant limits on the parameter t, if one looks 
only at leading terms. The coefficients of the leading terms are the 
peak values, and the only peak value (of holonomy, E , L, or angles) 
which depends on t is < L >. From eq. (jHSD even this peak value 
depends on t only via p/t, rather than p alone. 

<L + l/2>=p/t. (102) 

To put limits on t, one must consider the SC terms. I list various 
states contributing to the SC terms. First, the E SC terms, from 
eqs. (IHHj) and (|89|): 



v/l/(2 < L >) I ml(M)); v/l/(2t < L >2) | ml(L)). 

Each state is smaller than the leading term by the factor multiplying 
the state. Next, the holonomy SC terms, from eqs. (p9l) and (1 1001) : 

Vl/4 < L > I L + l,ml(M)); ^/tj2 \ L_) 
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One of these terms has t in the numerator, and one has t in the 
denominator, we can determine a best value of p and t by summing 
these two factors 

v/l/(2t < L >2) + v/t72, 

then minimizing the sum with respect to t. The resulting best value 
is 

t = 1/ < L > . (103) 

I cannot set t = l/< L > exactly, however, because then from 
eq. (11021) I must take p = 1. Appendix [Cl requires an expansion in 
exp(— p) ^ 1 to simplify D(H). As a compromise, I take p large, but 
not too large; say p = 5. Then the expansion of D(H) remains valid, 
since exp(— 5) is small; also the SC terms will be small, provided 
< L + 1/2 > is large enough. The t dependent factors suppressing 
the SC terms become 

v/l/(2t < L >2) = v/l/p<L>; ^/tj2 = ^J^l2 < L >. 

For < L > greater than 100 or so, and p = 5, these factors are 
sufficiently small. 

When the SC terms are taken into account, t is not arbitrarily 
adjustable. Values of p and t much different from p = 1 and t = 
l/< L > result in larger-than-optimal SC terms. 



E Fluctuations in Angles 

Given that the matrix h is strongly peaked at the value u it should be 
possible to translate this into a statement that the angles {9, 0) are 
strongly peaked at (a,/3), then obtain an estimate for the standard 
deviations of these angles. 

To begin, establish the two expectation values 



= {u,p\ (DW(/i) - D^'\u))oA I u,p)- (104) 

ordert,l/<L>= {u,p\ {D^^\h)] - D^^\u)])bo 

iD^'\h)-D('\u))oA\u,p). (105) 

The first equation is related to a mean; the second to a standard 
deviation. The first equation is zero because {D{h) — D{u)) \ u,p) 
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contains only SC states, and the leading SC states are odd in one 
of the variables, M or L— < L >, while {u,p \ is even. 

What happens when a second factor of D(h) - D(u) is applied 
to the state, as in eq. (11051) ? The right hand factor (D^^\h) — 
D^^\u))oA I UjP) is the sum of SC states given at eq. (1^ . 

\SCa) ■■= {D('\h)-D^'\u))oA\u,p) 

= N/N{L + l,ml{M))\L + l,ml{M)) 

[— cos fihA + sin fi{n x V)a] 
(2L + 1) 

-iN/N{L.) I L_)[sin/inA + cos/i(n x V)a]{106) 

The left hand factor {u,p \ {D^^\h)] — D^^\u)\) bo is the Hermitean 
conjugate {SCb \- The dot product of the two factors gives vector 
components of order unity times factors of order 

{N/[N(L + 1, ml(M))(2L + l)]}^ [N/N(L_)]2 

which are order 1/L and t, respectively, from eq. (llOOp .D 
We now have a sharpened version of eq. (11051) : 

{SCb I SCa) = {u,p\ p«(^)t - D^'\u)]bo) 

■{D^'\h)-D^'\u)),A\u,p), (107) 

with the left hand side given by eq. (11061) . 

The right hand side is a measure of the standard deviation of 
the holonomy. To turn this into a measure of angular standard 
deviation, expand D(h) around D(u). Define 



6 — a; 6(f) = (p — p. 



Then 



D{h)oA = d{a/2 + 66 /2)oAexp[tSzi(3 + 6(1)- Tr/2)] 

= d{a/2)oNexp[iSz{P - 'n-/2)]exp[i{69/2)S ■ h]NAGxp[iSz{6(t))] 
= D{u)oN[l + ^i6e/2)S -71 + tSzi6(l))UA- (108) 

The second line uses the behavior of the y axis under rotation around 
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exp[-iSz(/9 - 7r/2] SYexp[iSz(/3 - 7r/2)] = S ■ n. 
Now insert eq. (ITUg]) into eq. (ITOTD : 

(5Cb I SCa) = {u,p\ {D{u)[{69/2)S-h + Szm]}oA 

■{[{69/2)S-h + Szm]D{u)]}Bo I u,p). (109) 

All matrices D(u) have L = 1. On the right, the [66 6 (p) cross 
terms average to zero, leaving a positive definite expression linear 
in (< 59f >, < ((50)2 >. 

To separate the two standard deviations, first take A = B and 
sum over A. The matrix D(u) has its axis of rotation along n, there- 
fore 

D(u)S-hD(u)t = S-n. 

D(u) rotates the Z axis into the vector V, therefore 

D(u)SzD(u)t = S-V. 
We get 

l.h.s. = {N/[N{L + l,ml{M)){2L + l)]Y + [N/N{L^)]^- 
r.h.s = [<{6e/2y>+sm\a/2)<{6(j))^>]. (110) 

The left-hand side of eq. (I109P comes from eq. (I106P . The right-hand 
side uses 

[{S-Vr]oo = [L{L + 1)/2]{V^ + V^) 
= [L{L + I) / 2] sin^ {a/ 2), 

plus L = 1; and a similar formula for [S ■ n)"^, except n'j^ + Uy = 1- 
To separate (< 69^ > from < (50) ^ >), take A = B = in 
eq. (llOOp . This kills the 50 terms. In eq. (I106P the n terms disappear, 
since this vector has no Z component. The new left- and right-hand 
sides are 

l.h.s. = sm'^{a/2){N/[N{L + l,ml{M)){2L + l)]ysm'^H 

+ cos^ ft [N/N{L.)]^; 
r.h.s = [< {6e/2f > sin2(a/2)]. (Ill) 
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Eqs. flllOp and (11 111) may be solved for the individual standard 
deviations. 



< {59 /2f > = sin^ ;x(l/(8 < L + 1/2 >) 

+ cosV(V2)]; 
sm^{a/2) < (50)2 > = cos^ yu(l/(8 < L + 1/2 >) 

+ sinV(V2), (112) 

where eq. (llOOp was used to eliminate the ratios of normalizations. 
Eq. (11121) has the following geometrical interpretation. Sketch 
the following vectors: draw a vector V, in a direction {6/2, — 7r/2), 
very near the peak values (a/2,/5 — vr/2). Draw n, the axis of u, 
perpendicular to peak V and lying in the XY plane. Draw the 
third member of the orthogonal triad, n x V. Draw the angular 
momentum < L >. From the electron analogy, this is somewhere in 
the n, n X V plane, since it must be perpendicular to V. The exact 
location of < L > is tuned by adjusting the angle fi contained in 
the vector p: 

< E >~< L >=< L > (ua cos/i — h X V)a sin/x), 

from eq. Qj. 

Now suppose, for example, sin/x = 0, so that < L > is along n. 
The (t/2) term in eq. (11121) comes from fluctuations in magnitude 
L with the direction L fixed, while the l/<L + l/2> term comes 
from fluctuations in Lz (fluctuations in direction) with magnitude 
L fixed. (Cf. eq. (I85p : the standard deviation for fluctuations in 
L is yl/t), while the standard deviation for fluctuations in M is 
yL + 1/2.) Suppose first only L fluctuates, with the direction of 
< L > remaining fixed along n. Then from the sketch, as L—*L ± 
AL, the vector V will experience rotation along the direction 6: 
9 -^ 9±69. This checks with eq. (I112p : for sin/i = 0, a t/2 fluctuation 
(a fluctuation in the magnitude of L only) causes a fluctuation in 9 
only. 

Similarly, a 1/ < L+1/2 > fluctuation (a fluctuation in direction 
Lz, but not in magnitude L) causes a small fluctuation of L in a 
vertical direction. Since L must remain perpendicular to V, the 
latter fluctuates in an azimuthal direction. Eq. (11120 predicts an 
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azimuthal fluctuation, consistent with this picture. Other values of 
/i may be handled similarly. 

F SC States and Standard Deviations 

The uncertainty principle does not allow the SC terms to vanish. 
Consider, for example, the SC terms generated by the E operator. 
If they vanish, then the coherent state is an exact eigenstate of the 
spin. (The E operator brings down a factor of spin, therefore is 
essentially the spin operator.) If spin is exact, then from the un- 
certainty principle the canonically conjugate variables, the angles 
^,0), must be completely uncertain. This implies a completely un- 
certain holonomy, whereas a coherent state must have both spin and 
holonomy peaked. Therefore the SC terms cannot vanish. 

One can make a more quantitative statement about the SC terms 
and the uncertainty principle: the SC terms yield the standard de- 
viations of the operators E and /i*^^/^) around their means. Proof: let 

be an operator which is an approximate eigenvalue of the coherent 
state, 

Oa\u,p) = <Oa>\u,p) + SC 

= <Oa>\u,p) + ^eiA\i\nci). (113) 

i 

1 have given O a vector subscript because the important operators 
in this paper are vectors; but this feature is unimportant. The Cj 
are unit vectors. When the standard deviation is computed, the 
original coherent state drops out, and only the near coherent, SC 
terms contribute: 

{u,p\ {Oa- <OA>f \u,p) = y^^ej*^eiAX*Xi{nCj \ ncj) 

ij 

= ^Ci^ei^A-Ai, (114) 

i 

On the last line I have used the orthogonality of the near coherent 
states. From eq. (11141) . only the near coherent states contribute to 
the total standard deviation squared. 
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From the discussion of SC states in appendix [Dl e. g. eq. (IHHI) . 
each Aj will contain a factor of N/N(i). The above discussion clari- 
fies why each N/N(i) contains a factor which may be interpreted as 
a standard deviation, for either the L or M distribution. Eq. flll4p 
expresses the total squared standard deviation as the sum of indi- 
vidual squared standard deviations | Aj p contributed by each near 
coherent state. In statistics, the sum-of-squares form is obtained 
when the fluctuations in a total quantity are caused by fluctuations 
in several variables which are statistically independent. 
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